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Introduction 

An infrequent but significant hazard to life and property due to impacting asteroids and comets exists. 
Earth-approaching asteroids and comets also represent a significant resource for commercial exploitation, 
space exploration, and scientific research. The impact problem and those planetary bodies that could be a 
threat have been discussed in great depth in a wide range of publications. A popular planetary defense 
method is the deflection of asteroids and comets on a collision orbit with the Earth by changing their 
orbital velocities. It is fundamental to estimate required changes in the orbital velocities of dangerous 
objects to avoid a collision. 

A detailed optimization problem is formulated to calculate optimal impulses for deflecting Earth- 
crossing asteroids/comets, using nonlinear programming. The constrained optimization problem is based 
on a three-dimensional patched conic method to include the Earth’s gravitational effects and the asteroid/ 
comet’s orbital inclination. The magnitudes and impulse angles of optimal impulse change in velocity 
(AV) are accurately computed at various points on the asteroid/comet’s orbit to provide a given target 
separation distance. Based on these change in velocity analyses, we can establish the approximate cost 
and build strategies to prevent possible catastrophe due to the objects. 

Background 

The optimal AV and deflection strategy are dependent on the size and the orbital elements of the 
asteroid/comet, as well as the amount of warning time. Assuming a linear approximation between orbital 
energy and velocity increment, a velocity increment of about 1 cm/s is required to deflect an asteroid by a 
distance equal to 1 Earth radius for the order of a decade ahead of an impact (ref. 1). Using Keplerian 
motion and perpendicular impulse direction, both kinetic-energy deflection and nuclear-explosive deflec- 
tion are treated in reference 2. A method in reference 3 shows the instantaneous correction of asteroid 
velocity as a result of a spacecraft’s collision with an asteroid at perihelion only. In reference 4 the use of 
a nuclear explosive is discussed and it is concluded that a nuclear device having a yield of about 1 Mt 
may be required to deflect an object with a 0.3-km radius at a distance of 1 astronomical unit (au). Direct 
spacecraft impacts on Earth-crossing objects (ECOs) can provide kinetic energy to suitably disrupt up to 
about 0.1 -km stony asteroids and 0.3-km ice comets (ref. 5). One clear conclusion from these simplified 
analyses is that early detection gives longer reaction time, and asteroid/comet interceptions far from Earth 
are much easier and more desirable than interceptions near Earth because small deflections far away will 
produce a greater miss distance at the Earth. It is important to consider interception several orbital 
periods in advance as well as interceptions close to the Earth because both of these situations could be 
encountered depending on the size and orbital characteristics of an asteroid/comet. 

Recently, the astrodynamical optimization problem was formulated for Earth-crossing asteroids 
(ECAs) and presented accurate impulse angles as well as impulse magnitude according to impulse time 
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based on a two-dimensional, two-body analysis (ref. 6). The minimum impulse change in velocity (AV) 
requirement is not a monotonically decreasing function of warning time; rather, there is a finer structure 
associated with the orbital period of the colliding asteroid. Thus, the “optimal time” for application of AV 
is the earliest possible perihelion for warning times greater than 1 orbital period (ref. 6). However, this 
research has simplified the problem by assuming two-body orbital mechanics between the Sun and an 
ECA. This assumption neglects perturbations due to the Earth’s gravity. While these perturbations may 
not be present until the terminal phase of the impact scenario, they affect both long and short warning 
time analyses. Further research refines the heliocentric two-body analysis presented in reference 6 by 
including the gravitational effect of the Earth (ref. 7). The minimum AV is increased due to gravitational 
effects of the Earth, and effects on the minimum AV are sensitive to orbital elements of a target asteroid 
because the impact parameter is dependent upon those orbital elements. Gravitational effects of Earth 
have the strongest influence on minimum AV calculations for deflecting ECAs in nearly circular helio- 
centric orbits around the Earth. Hence, reference 7 concludes that the problem of deflecting ECAs should 
include the gravitational effects of Earth; the results in reference 7 followed from two-dimensional analy- 
sis. A near-optimal analysis for three-dimensional deflecting problems is described by using a system 
transition matrix (ref. 8). The study in reference 8 does not include Earth’s gravitational effects and is a 
first-order approximation because of the nature of the system transition matrix used. To have insights in 
deflecting Earth-approaching asteroids and comets with nonzero inclination, this paper expands the works 
of reference 7 by formulating an optimal deflection problem in three-dimension that also includes the 
Earth’s gravitational effects. The analysis of the impact deflection problem in this paper is based on a 
three-dimensional patched conic method. The analysis centers on how impulses applied to an asteroid at 
various points on the asteroid’s orbit affect the outcome when there is a presumption of collision other- 
wise. The analysis tool presented can be utilized in determining an accurate estimate for optimizing the 
time and position of intercepting the asteroid and comet for impact deflection. 

Deflection or orbit modification methods will be dependent on the size and composition of asteroids/ 
comets and the amount of warning time available. There are several deflection schemes such as the 
propulsive mode, kinetic deflection, nuclear detonation, and laser ablation. A study of spacecraft vertical 
landing on an asteroid indicates the possibility of giving a significant boost to an asteroid by means of a 
buried explosion of practical yield (refs. 9 and 10) and a high-thrust or low-thrust engine. Because every 
kilogram of propellant landed on an asteroid requires thousands of times more mass to lift it from the 
Earth and deliver it to the object, landing propellant on the asteroid in order to provide working mass is 
extremely inefficient (ref. 11). However, this strategy can easily control motion of an asteroid/comet 
because it has many opportunities to choose the direction and time of impulse application (ref. 3). Kinetic 
deflection and nuclear detonation demand much less “cost” than the propulsive mode and, consequently, 
require much less initial mass for deflection missions (ref. 12). A high-thrust, laser ablation technique 
should be seriously considered for use in deflecting or disrupting threatening celestial objects because 
these schemes appear both practical and possible. In this paper, we briefly discuss the potential of each 
strategy to deflect ECOs in conjunction with several future spacecraft concepts. 

Three-Dimensional Optimality Problem for AV Analysis 

Given an ECO and an established Earth collision, the problem is to minimize the impulse change in 
velocity (AV) required to deflect the ECO in such a manner as to miss the Earth by a minimum target 
miss distance (see fig. 1). The performance index is defined by the magnitude of the required deflection 
velocity: 


J(u) = 


AVj? + AVj5 + AV^ 


( 1 ) 
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TNW coordinate system 



Figure 1. Geometry of Earth-crossing asteroid/comet (not to scale) with TNW coordinate system. 


where AV T is the velocity increment aligned with the ECO velocity, AV N is the velocity increment 
normal to AV T in the ECO orbital plane, and AV W is the velocity increment normal to the ECO orbital 
plane. The impulse controls, u, are AV T , AV N , and AV W , which effectively give the magnitude 

and angle of the deflection velocity. The geometry of an ECA or comet 

with the “TNW” coordinate system is shown in figure 1. The calculation is always done to move the 
ECO’s trajectory from crossing the Earth’s orbit at the Earth’s center. This problem is subject to the 
heliocentric two-body equation outside the Earth’s sphere-of-influence (SOI) and geocentric two-body 
equations inside the SOI. Thus, the analysis of the impact deflection problem is based on a three- 
dimensional patched conic method. It is assumed that an ECO is influenced by the gravitational field of 
the Earth only when it is within the Earth’s SOI. Beyond the SOI, the object is considered to be affected 
only by the Sun’s gravitation. In the case of an impact scenario, the ECO begins in an elliptical orbit 
about the Sun. Once within the Earth’s SOI, the object’s motion is described by two-body orbit equa- 
tions for a hyperbolic orbit about the Earth. The radius of the Earth’s SOI is about 9.31 x 10 5 km, or 
0.00621 au. In figure 2, the approaching distance and a hyperbolic trajectory (ref. 13) are illustrated, and 
the geometry is still valuable for three-dimensional analysis. Under this framework, the constraints for 
the optimization problem can be described in terms of the terminal boundary conditions at the time (t f ) 
when the ECO intersects the Earth’s SOI: 


AV = 


AV^ + AV^+AV^ 


R-R SO i=0 (2) 

b-bi = 0 (3) 

R < 0 (4) 

where R is the distance between Earth and the ECO, Rsoi is the radius of Earth’s SOI, R is the time 
derivative of R, b is the approach distance of the ECO, and b ] is an impact parameter of the Earth. The 
constraints above mean that the approach distance ( b ) of the approaching ECO must be equal to the 
impact parameter (Zq) at the edge of the Earth’s SOI in order to deflect the ECO by a minimum target 
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Figure 2. Approach distance, b, at SOI of Earth (not to scale). 


miss distance R m . If the approach distance is less than the impact parameter, the ECO can be inside the 
target miss distance. When we set R m = 1 Earth radius, and b < b- v the ECO will collide with Earth. Note 
that if b = b v there will be a surface graze and the case is not considered as a collision in our analysis. If 
the approach distance is greater than the impact parameter, the impulse will not be a minimum value, 
forcing the ECO to miss the Earth by R m . The variable R is given by 


where r© is the radius vector from the Sun to the Earth, r a is the radius vector from the Sun to the ECO, 
\ a is the velocity of the ECO with respect to the Sun, and V 0 is the velocity of the Earth with respect to 
the Sun. In the heliocentric system, V„ is the vector difference at R S0I between V a and V© , so that 


The vector is described by its magnitude and an elevation angle, (p. From figure 2, the elevation 
angle is given by 



( 5 ) 


Then, we have 



( 6 ) 


V^V^-V© 


( 7 ) 


cos ((p + 90°) - Rsoj/V^Rsoi 


( 8 ) 


and the approach distance is obtained from 


Z7=R SO iCos (p 


( 9 ) 


The impact parameter, b \ , yields the expression 



(10) 
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where 


V esc = v/ 2 M-®/ R m (11) 

V esc is the escape velocity at R m , and U0 is the gravitational constant of the Earth. 

As noted in reference 6, closed- form solutions cannot be obtained for this optimization problem. To 
solve the heliocentric three-dimensional two-body motion, the method of Lagrangian coefficients (ref. 14) 
is used for describing r(t) and V(t) in terms of initial position vector, r(t 0 ), and velocity vector, V(t 0 ), in 
the inertial coordinate frame. At impulse time, t impulse , the original orbit of an ECO is perturbed by an 
impulse; hence, 


r (timpulse) r o(t impulse) 0 2 ) 

VOimpulse) — ^0 i m p U lse > AVxyz(t impulse) (13) 

AV X yz can he expressed in the inertial coordinate system (X,Y,Z) as follows: 

AVxyz = R 3 (-O) R 1 (-;) R 3 (— CO) R 3 (-a) AV TNW (14) 

where R 3 (-Q), Rj(-i), R 3 (-co), and R 3 (-oc) are rotation matrices about the longitude of ascending 
node (Q), the inclination (/), the argument of periapsis (co) for ECO, and the angle (a) between ECO 
perifocal coordinate and TNW system, respectively. Using r(t impulse ) and V(t impulse ) as initial 
conditions for the perturbed orbit, r(t f ) and V(t f ) are computed at the time when R = R S0I , If- 

Nonlinear programming (NLP) problems require determining x e R n , with x as the state vector, that 
minimize the scalar objective function J(u) subject to the equality constraints, Cj(x) = 0, and the inequal- 
ity constraints, Cj(x) < 0. It is assumed that the objective function is continuously differentiable through 
second order. The parameters needed to obtain the constraint equations can be calculated using the 
relationship between the orbital elements and the position and velocity vectors. For this deflection prob- 
lem, the free parameters are AV T , AV N , AV W , and t f , equations (2) and (3) are the equality constraints, Cj 
and C 2 , and equation (4) is the inequality constraint, C 3 . The parameter t f describes the time at which the 
constraints given by equations (2) through (4) are satisfied. The problem of minimizing the AV required 
for deflecting ECOs can now be cast in terms of a standard NLP problem. This formulation is also appli- 
cable to any Earth-approaching objects. 

Numerical Results of AV Analysis 

This analysis indicates the optimal impulses applied on an ECO at various points on the object’s orbit 
to provide the miss distance equal to R m when there is otherwise a presumption of collision at Earth’s 
center. The NLP formulated in the previous section was solved by using the MATLAB™ optimization 
tool box (ref. 15). In the discussions to follow, we will use the term impulse time to specifically mean 
either (tj mpact - ti mpu i se ) or its absolute value. Here ti mpact denotes the time at collision. Although the 
impact time is quite close to t f , it is not the same because t f is the time when R = R S0I . Also, because 
timpact an d ti m p U ise are not independent quantities, we choose ti mp act = 0. This initialization has the 
advantage of interpreting ti mp ulse as the time interval prior to impact if no action (i.e., AV maneuver) is 
undertaken. Also, it is apparent that we must have a warning time (i.e., the time interval between detec- 
tion and collision) greater than the impulse time. The impulse time provides a crude measure of the 
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warning time. The gravitational effects of the Earth are considered by using the three-dimensional optimi- 
zation problem described previously. The minimum required impulses for the deflecting problem have a 
targeted distance of 1 Earth radius. The solutions to the deflection problems represent impulse vectors 
that can be described by the magnitude of the minimum impulse and the optimal impulse angle. The 
impulse angle is described in the ECO’s orbital plane and is defined as the angle from the ECO’s original 
velocity vector to the impulse vector toward the Sun-ECO line, r a (see fig. 1). As reference 8 indicates, 
most numerical simulations show that the velocity increment (AV W ) normal to the ECO’s orbital plane is 
insignificant and ignored for all cases in these analyses. Thus, the magnitude of impulse consists of only 
AV t and AV N . For any given impulse time, the problem has two solutions for the optimal impulse angle 
separated by 180° while keeping the same magnitude of impulse (ref. 6). In this section, only one solu- 
tion for the optimal angle will be mentioned. 

There are two categories of impact scenarios: one (impact before perihelion, or preperihelial impact) is 
that an impact occurs before the ECO sweeps its perihelion (fig. 1), the other (impact after perihelion, or 
postperihelial impact) is that an impact occurs after the ECO passes its perihelion. The minimum impulse 
requirements for the two impact scenarios are similar but slightly different. The differences are caused by 
the varied geometric positions of the impulse point with respect to the Sun. Each impact scenario has two 
subscenarios: an impact occurs either at the ascending node or at the descending node of an ECO. These 
two subscenarios have identical impulse angles and magnitudes because the subscenarios have the same 
geometric positions of the impulse point. We present results from some AV analyses while our method is 
applicable to any ECO deflection mission. 

Near-Earth Asteroids 

There are three classes of near-Earth asteroids (NEAs): Atens, Apollos, and Amors. Aten-type aster- 
oids have a semimajor axis smaller than 1 au and an aphelion greater than 0.983 au, while Apollo-type 
asteroids have a semimajor axis greater than 1 au and a perihelion smaller than 1.017 au. Hence, Apollo- 
type and Aten-type asteroids can have Earth-crossing orbits. Amors have orbits that lie completely 
outside Earth’s orbit (perihelial distance between 1.017 and 1.3 au) but have the potential to be perturbed 
into Earth-crossing trajectories. 

We consider fictitious Apollo-type asteroids with semimajor axis a = 1.5 au, eccentricity e = 0.5, 
inclination i = 0°, 20°, 40°, and 60° first. For this example, the asteroids have orbital periods of approxi- 
mately 1.84 years, 0.75 au perihelial distance, and 2.25 au aphelial distance. Figures 3 and 4 include the 
magnitudes and angles of the optimal impulses. The impulse time is normalized to the period of the 
unperturbed asteroid for ease of interpretation. The abscissa represents the time when impulse is applied 
(as a fraction of the period of the asteroid) prior to collision. It can also be noticed that the separation 
between Earth and an asteroid, which can be achieved by an impulse, depends strongly on location of the 
impulse on the orbit as well as direction of the impulse with respect to the orbital velocity. The required 
minimum AV has a cyclic component imposed upon a secular variation that varies inversely with the 
impulse time. Generally, the minimum AV for inclined orbits (/ > 0°) is slightly less than that for planar 
orbits (/ = 0°). The case of preperihelial impact has more fluctuated magnitude variations than the case of 
postperihelial impact. Figure 4 contains a parallel history of the optimal impulse angle with respect to 
impulse time. It is apparent that the optimal angles for both planar and inclined cases are almost the same 
if the impulse takes place at more than half an orbital period of the asteroid before impact. There are 
fluctuations in the optimal impulse angle when the AV occurs less than half an orbital period before 
impact. The fluctuations are more for an orbit with a higher inclination than for one with a lower inclina- 
tion. Generally, postperihelial cases have more fluctuations in impulse angle than preperihelial cases. 
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Figure 3. Minimum AV: Apollo-type asteroids with a = 1.5 au, e = 0.5, i = 0°, 20°, 40°, 60 



Figure 4. Impulse angle: Apollo-type asteroids with a = 1.5 au, e = 0.5, i = 0°, 20°, 40°, 60 
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Figures 5 and 6 include the magnitudes and angles of the optimal impulses for another set of Apollo- 
type asteroids whose orbital elements are a = 2.0 au (orbital period =2.83 years), i = 20°, and e = 0.6, 0.7, 

0. 8, and 0.9. From these analyses, we know that higher eccentricity yields more fluctuated histories of 
impulse magnitude and angle. Figures 7 and 8 include the magnitudes and angles of the optimal impulses 
for a third set of Apollo-type asteroids whose orbital elements are e = 0.7, i = 10°, and a = 1.2 au (orbital 
period =1.31 years), 1.6 au (orbital period =2.02 years), 2.0 au (orbital period =2.83 years), and 2.4 au 
(orbital period =3.72 years). The minimum AV decreases with increasing orbital period (i.e., larger 
semimajor axis) because a longer orbital period yields a longer warning time. 

Figures 9 through 14 contain data demonstrating that the same trends are apparent for Aten-type 
asteroids. In figures 9 and 10, the semimajor axis and the eccentricity are fixed at 0.9 au (orbital period 
=0.85 years) and 0.4 au, respectively, and the inclination is varied from 0° to 60°. In figures 11 and 12, 
the semimajor axis and the inclination are fixed at 0.8 au (orbital period =0.72 years) and 10°, respec- 
tively, and the eccentricity is varied from 0.3 to 0.6. In figures 13 and 14, the eccentricity and the incli- 
nation are fixed at 0.6 au and 20°, respectively, and the semimajor axis is varied from 0.65 au (orbital 
period =0.52 years) to 0.95 au (orbital period =0.93 years). In all cases examined, the orbital inclination, 

1, has a relatively small effect on the minimum AV compared with eccentricity e and semimajor axis a 
(see figs. 7 and 1 1). 

When only a two-body approximation is applied to the deflecting problem, the minimum AV is line- 
arly proportional to the miss distance, R m (ref. 6). For example, in the two-body approximation, deflect- 
ing a dangerous celestial body by 10 R.0 requires exactly 10 times more AV than that needed to deflect 
the object by 1 R.0. When the gravitational effects of Earth are considered, the minimum AV is linearly 
proportional not to the miss distance but to the impact parameter, b r Figure 15 includes the minimum AV 
(impact after perihelion) of a fictitious asteroid whose orbital elements are a = 1.5 au, e = 0.5, and 
i = 20° as the miss distance R m varies from 1 R0 to 100 R s . One Earth-Moon distance is approximately 
60.27 R0. Analysis of equations (10) and (11) explains that the AV requirement for a deflection of NR0 
is less than N times the AV requirement for 1 R0 deflection. The reason is that V esc is reduced as R m is 
increased; hence, b i is not linearly proportional to R m . For instance, at t = -1 impulse time of the aster- 
oid’s orbital period, the AV requirement for 1 R0 deflection is about 2.5187 cm/s while the AV require- 
ments for 10 R @ and 100 R s deflection are about 21.6146 cm/s and 212.1871 cm/s, respectively. 

Short-Period Comets and Long-Period Comets 

Earth-crossing comets are classified into two types: short-period comets (SPCs), or those with an 
orbital period <200 years, and long-period comets (LPCs), or those defined here as having an orbital 
period >200 years. In figure 16, the histories of the minimum AV (miss distance of Earth radius) for a 
fictitious SPC whose orbital elements are given by a = 4.0 au, e = 0.85, and i = 0°, 20°, 40°, and 60° are 
shown. For this example, the comets have an orbital period of approximately 8.0 years, perihelial dis- 
tance of 0.6 au, and aphelial distance of 7.4 au. It is again noted that the AV W component is small enough 
to be neglected. The effect of inclination on the SPC minimum AV is also relatively small, as was 
observed previously for NEAs. Figure 17 includes the histories of optimal impulse angles corresponding 
to AV in figure 16. As impulse time approaches impact time, the class of postperihelial impact has peaks 
in optimal impulse angle while the class of preperihelial impact has an optimal impulse angle more sensi- 
tive to the inclination. Like other ECOs, the SPC has a slightly larger magnitude of minimum impulse for 
the case of preperihelial impact than for the case of postperihelial impact. For the SPC, when the impulse 
time is very close to impact time, the case of preperihelial impact also has more variation in the optimal 
impulse angle than the case of postperihelial impact as the inclination varies. The data in figures 18 and 
19 demonstrate that the same trends are apparent for LPCs, and the same trends may be expected for all 
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Figure 5. Minimum AV: Apollo-type asteroids with a = 2.0 au, i = 20°, e = 0.6, 0.7, 0.8, 0.9. 



Figure 6. Impulse angle: Apollo-type asteroids with a = 2.0 au, i = 20°, e = 0.6, 0.7, 0.8, 0.9. 
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a Preperihelial impact: local minimum occurred just before a minus integer impulse time 
b Postperihelial impact: local minimum occurred just after a minus integer impulse time 


Figure 7. Minimum AY: Apollo-type asteroids with e = 0.7, i = 10°, a = 1.2, 1.6, 2.0, 2.4 au. 
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Figure 9. Minimum AV: Aten-type asteroids with a = 0.9 au, e = 0.4, i = 0°, 20°, 40°, 60°. 



Figure 10. Impulse angle: Aten-type asteroids with a = 0.9 au, e = 0.4, i = 0°, 20°, 40°, 60°. 
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a Preperihelial impact: local minimum occurred just before a minus integer impulse time 
b Postperihelial impact: local minimum occurred just after a minus integer impulse time 

Figure 11. Minimum AY: Aten-type asteroids with a = 0.8 au, i = 10°, e = 0.3, 0.4, 0.5, 0.6. 
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a Preperihelial impact: local minimum occurred just before a minus integer impulse time 
b Postperihelial impact: local minimum occurred just after a minus integer impulse time 


. Minimum AV: Aten-type asteroids with e = 0.6, i = 20°, a = 0.65, 0.75, 0.85, 0.95 au. 



Figure 14. Impulse angle: Aten-type asteroids with e = 0.6, i = 20°, a = 0.65, 0.75, 0.85, 0.95 au. 
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Minimum AV, cm/s Minimum AV, cm/s 



Figure 15. Minimum AV: with respect to miss distance (R^). 



Figure 16. Minimum AV: short-period comets with a = 4.0 au, e = 0.85, i = 0°, 20°, 40°, 60°. 
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Figure 17. Impulse angle: short-period comets with a = 4.0 au, e = 0.85, i = 0°, 20°, 40°, 60°. 



Figure 18. Minimum AV: long-period comets with a = 40 au, e = 0.985, i = 0°, 20°, 40°, 60°. 
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Figure 19. Impulse angle: long-period comets with a = 40 au, e = 0.985, i = 0°, 20°, 40°, 60°. 


LPCs. In the figures, the semimajor axis and the eccentricity are fixed at 40 au (orbital period =253 years, 
perihelial distance =0.6 au, aphelial distance =79.4 au) and 0.985, respectively, and the inclination is 
varied from 0° to 60°. 

For ECOs, the data in figures 20 and 21 show AV and impulse angle histories of less than a 10-year 
impulse time. The AVs for the fictitious LPC (0.05 orbital period =12.7 years) are monotonically 
decreasing functions of impulse time, while the AVs for the fictitious NEA have the cyclic component 
varying with respect to orbital period. For a given miss distance and impulse time, the slower ECOs with 
respect to the Earth usually require a larger AV. The order of magnitude of AV is not significantly differ- 
ent for either the asteroids or the comets considered here. For an impulse time less than 10 years, the 
impulse angle of the LPC is mostly monotonically decreasing with respect to impulse time while the 
angle of the NEA has fluctuations due to its short orbital period. At the final stage, where impulse time 
approaches impact time, asteroids and comets have no significant differences in optimal AV and impulse 
angles. Additionally, the AV required for deflecting Earth-approaching asteroids and comets is dramati- 
cally increased, and impulse angle is quickly decreased. 

Interceptor Strategies 

It is important to estimate the required interceptor mass or energy to deflect or disrupt asteroids/ 
comets on a collision course with Earth. With the estimated mass or energy, we can investigate the 
appropriate strategies for the deflection missions. The interceptor parameter of M ht , Ei aser , M k , M nso , 
and M nss required for a given ECO mass and AV are expressed for high-thrust engine, laser ablation, 
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Figure 21. For comparison: impulse angle of postperihelial impact for various ECOs. 
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kinetic-energy deflection, standoff nuclear detonation, and slightly subsurface nuclear burst (refs. 3, 16, 
and 17) as follows: 


m AV 

M ht {kg} 


J laser 


m AV 

c m 


{J} 


. . m 

M t =- 


v e AV 

V v 2 


{kg} 


M, 


— mAV {kg} 
£{|) 


.. m v e AV 
A^nss ~ 0 , {kg} 

2e 4> 


(15) 

(16) 

(17) 

(18) 
(19) 


The interceptor parameter indicates mass (kg) of fuel for a high-thrust engine, energy (J) for a laser 
ablation technique, spacecraft total mass (kg) for a kinetic-energy deflection, mass (kg) of nuclear 
explosive for a standoff nuclear detonation, and a subsurface nuclear detonation. In these equations, other 
parameter values are defined or assumed as follows (ref. 17). The velocity of the ejecta is assumed as 

4 f dV 

— , density is 

v 2 ) 


v e = 100 m/s for subsurface bursts, m is the mass of the asteroid/comet, 


m = Ttp 

3 


assumed as 3 x 10 3 kg/m 3 for stony asteroids (200 kg/m 3 for comets), D is diameter of the asteroid/comet, 
about a fraction £ = 0.5 of the total energy could go into the kinetic energy of ejected material, nuclear 


explosions provide a specific energy <}> = 8 x 10 12 J/kg, a compression wave velocity in the asteroid/comet 
material is c p = 2 km/s, c e is the exhaust velocity of a high-thrust engine, c m is the laser momentum 
coupling coefficient (assumed in ref. 1 6 as c m = 5 dyn-s/J), v is the orbital velocity of an asteroid/comet, 
and AV is the orbital velocity increment to deflect an asteroid/comet on a collision course with Earth by 1 

Earth radius. It is assumed that the velocity of the interceptor is much less than that of celestial objects 
and can be neglected. 


The optimization method mentioned previously yields the AV to be used in equations (15) 
through (19). Note that equations (15) through (19) implicitly have estimation errors to give the final 
interceptor parameter, and the impulse is assumed to be applied with the optimal impulse angle discussed 
previously. Two fictitious asteroids (Aten-type asteroid: a = 0.9 au, e = 0.4, i = 0°, orbital period = 0.85 
years; Apollo-type asteroid: a= 1.5 au, e = 0.5, i = 0°, orbital period = 1.84 years) and two comets (SPC: 
a = 4.0 au, e = 0.85, i = 0°, orbital period = 8.0 years, LPC: a = 40 au, e = 0.985, i = 0°, orbital 
period = 253 years) are used. With AV calculated with respect to impulse time, equations (15) 
through (19) can be applied to estimate the final interceptor parameter to deflect the asteroid/comet by at 
least 1 Earth radius. Inclined asteroid/comet orbits (/' ^ 0°) should give approximately the same intercep- 
tor parameter as that of a planar orbit because there are no significant differences in AV histories for both 
the inclined and planar orbits. 


Figures 22 and 23 include the histories of the final interceptor parameter in order to deflect 1-km 
diameter stony asteroids (Aten-type or Apollo-type) by 1 Earth radius, while figures 24 and 25 include 
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Interceptor parameter, kg Interceptor parameter, kg 



Figure 22. Various interceptor masses for 1-km stony Aten-type asteroid. 



Figure 23. Various interceptor masses for 1-km stony Apollo-type asteroid. 
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Time before impact when mitigation is applied, years 
Figure 24. Various interceptor masses for 1-km short-period comet. 



Time before impact when mitigation is applied, years 


Figure 25. Various interceptor masses for 1-km long-period comet. 
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similar histories for 1-km diameter comets (an SPC and an LPC). The final intercepter parameter for 
Apollo-type asteroids shows the same pattern as those for Aten-type asteroids, with less fluctuations with 
respect to impulse time. For these figures, the payload masses of nuclear equipment are estimated in 
table 1, which aligns the data to show the relationship between nuclear yield and nuclear device weight 
(ref. 18). Note that a minimum payload mass of 50 kg is assumed when nuclear yield is less than 10 kt. 
The final interceptor parameters also vary according to the fluctuations of minimum AV required for the 
space deflection missions. For both the asteroid and comet, a high-thrust method requires much more 
mass than other methods identified (kinetic-energy deflection, standoff nuclear detonation, and subsurface 
nuclear burst) because of the fuel expenses for a soft landing and smaller exhaust velocity (ref. 3). 
Kinetic-energy deflection methods require hundreds of times less mass than high-thrust methods for 
asteroid deflection, and several orders of magnitude less mass for comet deflection. Nuclear detonation 
methods reduce the final payload mass by several orders of magnitude compared with the kinetic-energy 
deflection method. The explosive mass of a slightly subsurface nuclear burst is less than that for standoff 
nuclear detonation by several factors of 1 0 for both asteroid and comet deflection. However, the payload 
mass of both nuclear methods could be 50 kg when impulse time is greater than 1 year before impact. An 
interesting phenomenon is that, for the SPC case, spacecraft mass for kinetic-energy deflection could be 
less than nuclear payload mass (fig. 24) when deflection is provided at the AV’s local minimum. The 
reason for this phenomenon is that the velocity of the comet is relatively high at perihelion, and final mass 
of kinetic deflection is inversely proportional to the square of orbital velocity. Generally, the interceptor 
masses for comet deflection are less than those for asteroid deflection, because densities of comets are 
assumed to be much less than those of asteroids. Because the mass of an asteroid and comet is linearly 
proportional to the cube of its diameter, it is easy to estimate the final interceptor masses for different 
sizes of asteroids and comets. For example, 0.1 -km stony asteroids require one thousand times less final 
mass than 1-km stony asteroids. For asteroid missions, the “optimal time” for application of the AV is the 
earliest possible perihelion for impulse times greater than 1 orbital period, while for LPC missions, the 
optimal time is simply the earliest possible time. Surprisingly, for an LPC mission, the interceptor mass 
for the kinetic-energy deflection is approximately the same when the impulse time goes from 1 to 
10 years before collision with Earth, although AV is increased during the period (see fig. 18). The reason 
for this phenomenon is that the orbital velocity of the comet is increased accordingly as it travels toward 
its perihelion, and the interceptor mass of kinetic-energy deflection is inversely proportional to the square 
of orbital velocity. 


Table 1. Approximate Nuclear Yield Versus Nuclear Device Weight 3 


Nuclear yield, kt 

Payload weight, kg 

10 

50 

100 

100 

1000 

1000 

10000 

20000 


a From reference 18. 


Use of a Hybrid Propellant Module 

A hybrid propellant module (HPM) is a future reusable in-space transportation concept being studied 
under NASA’s Revolutionary Aerospace Systems Concepts (RASC) Program. The HPM can store 
indefinitely both chemical and electrical propellants and provide propulsion with attached modular orbital 
transfer/engine stages. The HPM can utilize a chemical propellant (liquid hydrogen/liquid oxygen) for 
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fast orbital transfers and can use an electrical propellant (low thrust) for prepositioning or to return the 
HPM for reuse and refueling. A chemical transfer module (CTM) serves as a high energy injection stage 
when attached to the HPM, while a solar electric propulsion (SEP) module serves as a low-thrust transfer 
stage when attached to the HPM. Although propulsive AV is much less efficient for deflecting asteroids 
and comets, the cost of utilizing an already available infrastructure might be favorable compared with 
developing a separate defensive capability. A preliminary investigation of how an HPM-based transpor- 
tation system could be utilized is provided. It is assumed that the dry mass of an HPM is 4000 kg, the dry 
mass of a CTM is 4400 kg, the chemical propellant mass of a CTM is 30000 kg, and the effective exhaust 
velocity, c e , is 4.4 km/s for the CTM. An HPM/CTM/SEP vehicle stack can travel to an ECO by using 
the low-thrust, high-specific impulse provided by the SEP and the high thrust provided by the CTM and 
HPM chemical propellant. The idea of having abundant propellant sources in deep space may offset the 
disadvantages of a chemical propellant when combined with electric propulsion. 

After rendezvous with the ECO, the combined vehicle could use its CTM as a high-thrust engine with 
a full 30000 kg load of chemical propellant to change the orbit of the asteroid/comet. According to 
figures 22 through 25, by using a CTM as a high-thrust engine, it can be estimated that =0. 1 3-km asteroid 
(=0.3-km comet) defense mission could be performed 2 years before collision. The same CTM could be 
used for =0.06-km stony asteroid (=0. 1 5-krn comet) deflection just 6 months before collision. An asteroid 
of about a 1-km (0.2-km) diameter can be deflected by direct spacecraft (HPM with CTM, total dry mass 
is about 8400 kg) impact on the target for longer than a 2-year (6-month) impulse time if the impulse is 
applied to an asteroid at its perihelion. An approximate 1-km comet can be deflected by kinetic energy 
using HPM with CTM for longer than 6 months of impulse time. These cases assume that all the fuel of a 
combined HPM with CTM can be consumed to travel longer or faster to dangerous objects. With an 
impulse time of longer than 2 years, HPM (dry mass of 4000 kg) itself could perturb orbits of asteroids/ 
comets of 0.3 km in size using direct kinetic impact on the objects. By assuming that the HPM/CTM 
stack could deliver a 20000-kg nuclear warhead, this in-space transportation system could be used for a 
1-km asteroid or a 2-km comet deflection mission for any short impulse time. 

An HPM-based in-space transportation system could be available in 10 to 20 years. For comparative 
purposes, it is useful to determine what deflection strategies are currently possible. Payload mass is 
directly related to propellant mass, so the final mass of an interceptor can be converted into a rough 
estimate of interceptor cost. An interceptor payload would require 3 times more mass into deep space and 
10 times more mass into low-Earth orbit (ref. 19). The Russian Energia launch vehicle was the most 
powerful launch vehicle and had a payload capacity of 80000 to 90000 kg for low-Earth orbits (ref. 12). 
With this launch capability, about 9000 kg of maximum final interceptor mass can be reasonably 
assumed. The data in figures 22 through 25 indicate that a 0.1 -km asteroid or 0.3-km comet defense 
mission could be performed for any short impulse time by using kinetic-energy technology. A 1-km 
asteroid would require at least 2 years of impulse time if the AV is applied at the asteroid’s perihelion, 
while a 1-km comet would require at least 1 year. These results coincide with previous studies stating 
that kinetic-energy impacts are suitable to deflect rocky asteroids of up to about 0. 1 km in size (ref. 5). It 
was also shown that today’s technology is sufficient to defend asteroids of less than 0.1 km in diameter 
for warning times on the order of several years (ref. 12). A 1-km asteroid or a 2-km comet defense mis- 
sion could be prepared for a few months impulse time by using nuclear detonations. Nuclear explosives 
offer energies great enough for much larger asteroids/comets and shorter impulse times. The studies in 
this paper may be slightly different from previous studies due to different approximations and assump- 
tions. Generally, for any warning time, we can roughly say that the kinetic-energy deflection would be 
preferable for smaller ECOs (less than a few hundred meters in size); the nuclear strategy would be 
required for large objects (greater than 1 km in size) or for very short impulse times. 
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Use of a Laser Ablation Technique 

Although nuclear detonations and kinetic-energy impacts appear to be the most practical near-term 
impact protection options, both methods have uncertainties in the momentum change imparted and can 
potentially fragment the asteroid/comet. This fragmentation could create multiple impacts, which ulti- 
mately could be more hazardous than the original object. Additionally, it is envisioned that comets and 
asteroids could provide a virtually limitless supply of resources for the future exploration and colonization 
of space. An orbit modification technique that could also facilitate the utilization of nonthreatening 
comets and asteroids would likely function in a controlled manner. One approach to alter the trajectory of 
the object in a highly controlled manner is to use pulsed laser ablative propulsion. This technique has 
been studied to remove space debris from low-Earth orbit (ref. 20). A sufficiently intense laser pulse 
ablates the surface of the space debris by causing plasma blowoff. The momentum change from a single 
laser pulse is very small. However, the cumulative effect is very effective because the laser can yield 1 0 
to 1000 pulses/s over several minutes (ref. 21). The dynamic reaction from multiple laser hits reduces the 
perigee altitude of the orbiting debris. Once perigee is reduced sufficiently (=200 km), atmospheric drag 
deorbits the debris. The highest pulse energy currently available is about 20 kJ/pulse operated at under 
0.02 Hz (ref. 20). In the future, a much more advanced version of this same technique may be used for 
the deflection/orbit modification of asteroids/comets on a collision trajectory with Earth. The laser abla- 
tion technique could overcome the mass penalties associated with other nondisruptive approaches because 
no propellant is required to generate the AV (the material of the celestial object is the propellant source). 
Additionally, laser ablation is effective against a wide range of surface materials and does not require any 
landing or physical attachment to the object. 

For deflection of distant asteroids and comets, the power and optics requirements of a laser ablation 
system on the ground or near Earth may be too extreme to contemplate in the next decades. One hybrid 
solution would be to permit a spacecraft to carry a laser as a payload to a particular celestial body. The 
spacecraft would require an advanced propulsion system capable of rapid rendezvous with the asteroid/ 
comet and an extremely powerful energy generation system. The spacecraft would orbit or station-keep 
with the object at a “small” standoff distance. In this section, it is assumed that a spacecraft with a laser 
ablation tool has already rendezvoused with the asteroid/comet. Figure 26 includes the required laser 
energy for various ECOs to deflect them by 1 Earth radius; it is demonstrated that the required laser 
energy also varies according to the required AV. The size of the celestial bodies in figure 26 is assumed 
to be 1 km. Again, it is easy to estimate the final interceptor energy for different sizes of ECO because 
the energy of an ECO is linearly proportional to the cubic of its diameter. Figure 27 houses an explana- 
tory detail of figure 26 for an Apollo-type asteroid and an LPC, using fixed impulse times, and describes 
the required energy for laser ablation to deflect the given size of an asteroid and comet by 1 Earth radius. 
Figure 28 includes the achieved energy for a given laser power and operation period. From figures 27 
and 28, we can easily estimate the power of laser for a selected case. For instance, if we can use a 1-year 
impulse time, about 2.5 x 10 3 GJ of energy (fig. 27) is required for a 0.1-km stony asteroid to be deflected 
by 1 Earth radius. The data in figure 28 can be used to determine what laser power is required and how 
long the laser should be continuously operated in order to achieve 2.5 x 10 3 GJ of energy. If we choose a 
laser system with 1 MW of power (100 kJ/pulse and 10 Hz laser repetition frequency), it would take 
about 30 days of operation to provide 2.5 x 10 3 GJ of energy. If we choose a laser system with 10 MW of 
power, it would take about 3 days of operation to provide 2.5 x 10 3 GJ of energy. If a laser system with 
average power of 1 MW could be continuously used for asteroid deflection 2 years before collision with 
Earth, it would take about 8 days to provide a 0.1-km stony Apollo-type asteroid with enough energy 
to avoid the collision. If the same equipment is available 6 months before the collision, it would take 
60 days for the same asteroid. For deflection of a 1 -km LPC, a 25-MW laser system needs to be operated 
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Figure 26. Required laser energy of postperihelial impact for 1-km ECOs. 



Figure 27. Required laser ablation energy versus size for Apollo-type asteroid and long-period comet. 
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Figure 28. Achieved energy for given laser power and operation period. 


continuously for about 3 months beginning 1 year before collision. Although these power levels may not 
be readily achievable today, they are consistent with the power levels needed for high-thrust, high- 
specific impulse propulsion systems envisioned for rapid in-space transportation in the future. This 
power needed for distant rendezvous with an asteroid/comet would then be available for the orbit modifi- 
cation phase of the defense mission. 

Overview of Results 

The minimum AV vector has finer structures associated with the impulse time to deflect an ECO. 
Both asteroids and comets have similar fluctuation structures of AV and impulse angle histories with 
respect to their orbital periods. For a given miss distance and impulse time, the slower Earth-crossing 
objects relative to the Earth usually require a larger AV. When the gravitational effects of Earth are 
considered, the minimum AV is linearly proportional not to the miss distance but the impact parameter, b v 
The effects of orbital inclination on AV and impulse angles are very small for the deflection problem. 
The AV and impulse angles are more sensitive to eccentricity and orbital period rather than the inclina- 
tion. The velocity increment, A V w , normal to orbital plane of Earth-crossing objects is negligible for 
most deflection missions. The case of preperihelial impact has slightly larger local minima in the magni- 
tude of impulse than the case of postperihelial impact, and preperihelial impact cases have more fluctua- 
tions in optimal impulse angle than postperihelial impact cases. At the final stage, where impulse time 
approaches impact time, the optimal impulse has been changed dramatically and the AV is logarithmically 
increased. Earth-crossing asteroids and comets may have no significant differences in optimal AV and 
impulse angles when impulse time is within a few months. 
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The final interceptor mass (or energy) also changes according to the fluctuations in the minimum AV 
required for the space deflection missions. The optimal deflection strategy is dependent on the size and 
composition of an asteroid/comet as well as the amount of warning time. The potential of several deflec- 
tion approaches, including high-thrust engine, kinetic deflection, nuclear detonation, and laser ablation, is 
discussed. A high-thrust method requires much more mass than other methods. Kinetic-energy deflec- 
tion needs about a few tens or a few hundreds times less mass than a high-thrust method. Nuclear deto- 
nations reduce final payload mass by several orders of magnitude compared with the kinetic-energy 
deflection method. For a given diameter target, the interceptor mass for comet deflection is usually less 
than the mass needed for asteroid deflection. Although not as mass or energy efficient as a nuclear deto- 
nation or kinetic-energy impact, a laser ablation system has significant advantages and the potential to 
provide controlled orbit modification of Earth-crossing celestial bodies. 

Concluding Remarks 

Deflection or orbit modification of potentially Earth-impacting asteroids/comets has been considered. 
Several techniques were assessed. 

These analyses demonstrate that the optimal deflection strategy is highly dependent on the size and the 
orbital elements of the asteroid/comet, as well as the amount of warning time. It is concluded that the 
best deflection strategy and optimal impulse change in velocity (AV) should be carefully investigated if 
an impacting object is discovered, provided that sufficient time is available. Depending on when the AV 
can be applied to an Earth-crossing object (ECO), the optimal impulse angle is not necessarily along the 
velocity or antivelocity vector, particularly when the impulse time is less than 1 orbit before impact. 
Although not as mass or energy efficient as a nuclear detonation or kinetic-energy impact, a laser ablation 
system has significant advantages including the potential to provide controlled orbit modification of 
ECOs. Regardless of the method employed, early application of the impulse dominates the magnitude of 
the deflection effort. Impacting objects confirmed years or decades in advance are much easier to deflect 
and allow a detailed assessment and development of a planetary defense strategy. However, the use of 
operational space-based elements in the event of a confirmed threat could provide a reliable, rapid 
response to ECOs, even with relatively short warning times. 
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